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Abstract. The decays of light vector mesons into three pseudoscalar mesons are calculated to leading order 
in the recently proposed counting scheme that is based on the hadrogenesis conjecture. Fully differential 
as well as integrated decay widths are presented. Since the required parameters have been fixed by other 
processes, the considered three-body decays are predictions of the presented approach. The decay width 
of the omega meson into three pions agrees very well with experiment. The partial decay widths of the K* 
into its three K-ktt channels are predicted. 



1 Introduction 

One of the open challenges of QCD is to develop a system- 
atic scheme for the calculation of hadronic reactions and 
decays. In the sector of mesons made out of light quarks, 
chiral perturbation theory rjj[2] provides such a system- 
atic approach. However, it is restricted to energies where 
no other mesons than the Goldstone bosons participate 
as active degrees of freedom. It is therefore important to 
construct effective field theories with additional degrees of 
freedom. In particular, in the framework of the hadrogene- 
sis conjecture [3] such an extension is mandatory since the 
pseudoscalar and vector mesons together with the baryon 
octet and decuplet are regarded as the preselected states 
which build up the other hadrons — at least in the low- 
energy regime. 

Recently we have suggested in [4] a counting scheme 
for flavor-SU(3) systems of Goldstone bosons and light 
vector mesons. The quest is to have a scheme at hand 
which allows in particular for decays of vector mesons into 
pseudoscalar states. The masses of pseudoscalar (P) and 
vector mesons (V) and all involved momenta are treated 
as soft, i.e. mp, my ~ Q, where Q is a typical momen- 
tum. Such an assumption may be justified in QCD with 
a large number of colors N c , which implies that the chiral 
symmetry breaking scale is much larger than the vector 
meson masses, i.e. 4^/^ 3> my. Consequently, all chirally 
covariant derivatives D^ are ~ Q, no matter whether they 
act on pseudoscalar or vector mesons. This is in contrast 
to the approach where vector mesons are treated as heavy 
matter fields [5]. The latter formalism does not allow for 
a systematic treatment of processes where the number of 
vector mesons changes, like e.g. in vector-meson decays. 
Our scheme is also different from the one presented in [6l[7] 
where all derivatives are treated as soft, but the masses of 
the vector mesons are treated as heavy. Strictly speaking, 
such an approach is restricted to energies below the vector 



meson masses. Nonetheless, resonance saturation models 
based on large-7V c considerations have been built upon 
the works [6l[7j and have been applied also to processes at 
higher energies, see e.g. [8] and references therein. 

For the present purpose it is only important to recall 
from [3] that loops are suppressed which is a generic con- 
sequence of large- N c considerations [9. It is important 
to note that the vector mesons are represented by tensor 
fields in our approach. The chosen representation cannot 
change the physical results, but contact terms (e.g. four- 
point interactions) can look different in different repre- 
sentations and different orders in a given counting scheme 
might be attributed to them. 

It has been demonstrated in ;4 a that the proposed 
scheme can be used to describe hadronic and electromag- 
netic two-body decays of vector mesons. In the present 
work we explore some further consequences of this frame- 
work by calculating to leading order the decays of vector 
mesons into three pseudoscalar mesons. It will turn out 
that all necessary parameters for the leading-order calcu- 
lation have already been determined from the two-body 
decays. Therefore, our calculations can be tested against 
the available experimental data and we can provide pre- 
dictions for processes which were not measured yet or not 
measured in a fully differential way. 

In the following we will study the process lo — > in 
and the three decay branches of K* — > Kirir. There are 
no data for any differential widths of these three-body 
decays available. The partial decay width for the omega 
decay is measured QJj] whereas for the K* decays only 
upper limits exist [IT]. The omega decay has been stud- 
ied by several groups, see e.g. [12] and references therein. 
Typically, however, the corresponding coupling constants 
have been fitted to the partial decay width of the process 
lo — > 3n. In our approach we have determined the required 
coupling constants from two-body processes and we have 
a counting scheme which tells us that we have covered 
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all relevant vertices. An approach which is on the tech- 
nical level close to ours is the one presented in [8]. Also 
there tensor fields have been used to represent the vector 
mesons. However, the counting scheme is different since 
vector meson masses are treated as heavy in [3]. Therefore, 
also the results cannot be compared easily. For the decays 
of the K * the literature is scarce. In [13] various models 
including vector-meson dominance and also contact inter- 
actions (see discussion below in sect. [2]) have been applied 
to the three-body decays of the K* . We will comment on 
the work of (13j when we present our results below. In con- 
trast to previous works differential decay widths for the 
considered processes will be presented. To the best of our 
knowledge this issue has not been studied at all by other 
theory groups. 

The work is structured in the following way: In the 
next section we present the relevant part of our leading- 
order Lagrangian and provide some general formulae im- 
portant for three-body decays. In sect. [3] we calculate the 
decay of the omega meson into three pions. We will present 
the fully differential decay width (Dalitz plot) as well as 
integrated quantities thereof. The corresponding calcula- 
tions for the decay channels of the K* into one kaon and 
two pions are presented in sect. [U Finally in sect. [5] we 
summarize our results. 



2 Three-body decays 

We compute the decays of light vector mesons into three 
pseudoscalar mesons. In general, such a calculation in- 
volves (a) tree diagrams with three-point vertices, (b) a 
tree diagram with a four-point vertex and (c) loops. We 
shall show now that only case (a) is relevant at leading 
order in the scheme proposed in [3]: For the considered 
decays the leading-order calculation is of second order. 
Loops are suppressed contributing at fourth order [4j. In 
principle, it is conceivable that a four-point vertex with 
one vector and three pseudoscalar states could enter at 
leading, i.e. second, order. However, within the tensor rep- 
resentation, which we use, one needs at least four deriva- 
tives to construct such a vertex. Thus it has at least or- 
der four according to the counting rules of [I]. Hence at 
the leading-order calculation, on which we concentrate in 
the following, one only needs three-point vertices. Since a 
three-point vertex with three pseudoscalar states does not 
exist, the decay process happens in two steps: The initial 
vector meson decays into another vector meson and one 
pseudoscalar state. The emerging vector meson finally de- 
cays into two pseudoscalar states. The generic diagram is 
depicted in fig. [TJ 

We recall from [4] the relevant part of the leading-order 
Lagrangian: 
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tr {V 1 "' dp$ d u $} 



Fig. 1. Feynman diagram which contributes in leading order to 
the decay of a vector meson (solid line) into three pseudoscalar 
mesons (dashed lines). 



with the conventions 
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e^ Q/3 tr{[V,V/ Q/3 ] + [^,xo] + }. (1) 



In ([T]) the quantity / denotes the pion-decay constant 
in the chiral limit. The quantity my has been introduced 
for convenience to render the coupling constants in |T]) di- 
mensionless. The combination my hp/ f 2 has been fixed in 
[4] from the decays of vector mesons into two pseudoscalar 
states. The leading-order Lagrangian of our approach has 
only one such term. Hence it predicts a universal value for 
the coupling constant for all these decays, i.e. for p — > 27r, 
K* — > ir K and <fi — > K K. Indeed, the values for hp ob- 
tained from fits to the experimental decay widths agree 
with each other within ±10%. The values for Ha and &a 
have been determined from the radiative decays of vector 
mesons [4]. Here, it turned out that one universal vertex 
would be insufficient to describe all the radiative two-body 
decays of the vector meson nonet. Indeed, our leading- 
order Lagrangian yields two terms — one flavor symmetric, 
one flavor breaking — which are then sufficient to describe 
the data [4] . We regard this as a support of our approach. 

The only decays which are allowed by energy-momen- 
tum conservation and where our Lagrangian yields non- 
vanishing couplings at leading order are ui — ► 3-7T and 
K* — > 2ir + K. Note that in leading order the <fi is purely 
strange and OZI suppression [T4l[i31[T6] is perfect. Hence 
we do not provide results for the OZI suppressed decay 
branch <j> — ► 37r. 

In the following we use the values 

/ = 90 MeV, m v = 776 MeV, 

h P = 0.304, h A = 2.l, b A = 0.27. (5) 

Note that the precise value for hp emerges from a fit to 
the electromagnetic pion form factor [17]. For the two- 
body decay of the K* this value for hp implies a partial 
width of 50.9 MeV, in excellent agreement with data [10] . 
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In [4] we have deduced for hp a value of 0.29 ± 0.03 which 
agrees with the one given in (|5|) within the error bars. 
This result of [4] emerged from a simultaneous fit to the 
decay widths p — > 7r n, K* — > K it and <j> — > K K. For the 
processes of interest in the present context the 4> decay 
does not enter. Hence we use a value for hp which agrees 
well with the processes p — > n n and K* — > K it relevant 
for the decay chains depicted in fig. [TJ 

We consider the decay of a state V with momentum p 
into three states with momenta p\ , P2 and p 3 . According 
to [TU] the double differential decay rate 



one parameter combination ~ hp hA determines the de- 
cay process, at least in the chiral limit. Quantitatively the 
extra term ~ m 2 - bA plays only a minor role. 

In fig. [2] the Dalitz plot for the decay process uo — ► 
7T+ 7r~ 7T° is depicted using the values (JSJ). The shape can 
be compared with the one obtained from pure phase space, 
i.e. replacing 0^^,^+ „.- ,« by a constant. A pure phase 
space Dalitz shape is shown in fig. [3J There are no visible 
differences and indeed in the kinematically allowed range 
the squared matrix element |C( i) _^ 7r + „■- w o| 2 is rather flat 
(not shown). 



V^123 



1 



dm\ 2 dm\ 3 (2tt) 3 32 M 3 



a 



V^12 3 



P 



■t^ap^ ti.titpsaPlPti.fi' (6) 



is determined by a matrix element Cv—>i 2 3 and the phase- 
space factor P. In © we introduced the following variables 



= (pi+Pjf 



M 2 = p 2 , 

4 = TO 12 , 



2 2 2 2 2,2,2,2 /»\ 

q 2 =m 13 =p -m 12 -m 23 +p 1 +p 2 +p 3 . (7) 

The additional quantities g 2 are introduced for later con- 
venience. 



3 Decay 



U! 



3tt 



The differential decay width of the process u> 
is determined by the matrix element 
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Fig. 2. Dalitz plot for the decay process lo 



Spi^l) ill + m l>) + Sp(ql) (q 2 + m^)] trary normalization). Note m\ 2 



0.5 



(arbi- 



and rrin 



2 my hp m\ h A 



[sM) 



where we use the kinematics of 
propagator 

S P (q 2 ) = 



SMI) 
0. In 



S P (ql) 



(8) 



q 2 — m 2 



^ p {q 2 ) 



the rho-meson 



(9) 



The single-differential decay width 



dr^ 



,3tt 



dm 2 2 



is plotted 



appears. Here r p (q 2 ) denotes the energy-dependent width 
of the rho meson. It turns out that the results only change 
by about 2%, if this width is neglected. We do not think 
that our leading-order calculation and the determination 
of our coupling constants have a comparable accuracy. 
Hence we refrain from a more detailed modeling of the 
rho-meson propagator including e.g. an energy dependent 
real part of the self energy. 



The factors qf 



which accompany the rho-meson 



propagators in ([5]) emerge from the use of the tensor repre- 
sentation for the vector mesons. In other representations 
this combination might split up into contact terms and 
vector-meson exchange terms. In our approach we do not 
have extra contact terms in leading order and therefore 



in fig. [4] One observes the typical rise and fall of phase 
space. We predict that the shape of the fully and single- 
differential decay widths do not differ significantly from 
pure phase space. It was checked that this finding emerges 
from the interference of the respective three rho-meson 
terms which appear in (jSJ) . Physically this corresponds to 
the three decays u> — ► p + ir~ , u> — > p°ir°, to — > p~ ir + . 
If there was only one rho meson, e.g. the p + , the Dalitz 
shape would differ significantly from pure phase space. 
For the partial decay width one gets 

r^ 37r = 7.3MeV (10) 

which agrees very well with the experimental value [10J 

r^ 3w = (7.57 ±0.13) MeV. (11) 

As stated before the influence of bA is of minor importance 
here. Using bA — leads to r u ->an — 7.7 MeV. We empha- 
size that our coupling constants have been obtained from 
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Dalitz plot for the phase space of the decay process 



which connects in a systematic way the two-body decays 
studied in 0] with the three-body decays studied here. 



4 Decay K* 



K Till 



We turn to the decay modes of the K* + . Due to isospin 
symmetry the decay modes of the K* a need not be con- 
sidered separately. Three decay channels exist: K* + — > 

7T+ 1t~ K+, K* + -> 7T° K+ 7T°, K* + -> 7T+ K° 7T°. 

Let us first discuss the experimental situation: In [11] 
upper limits are given for the processes 



r K -^ + n-K- <40keV, 
< 35keV, 



r 



K* 



*TV~ 7T° K° 



r R*«^+K-K° < 35keV. 



(12) 



Using particle-antiparticlc symmetry and flipping the iso- 
spin for the latter process this yields the following con- 
straints for the K* + decays: 



r K , + ^ +7T - K+ <40keV, 
r K *+^-K+T!°K° < 35keV, 
r K .+_^+^- K + < 35keV. 



(13) 
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Fig. 4. Single-differential decay width for the process u — » 
7r + 7T~ 7T°. The solid line gives the full result, the dotted with 
b A = 0. 



other processes. Thus, the good agreement between (jTOJ) 
and (jlip can be regarded as a further justification for our 
approach which predicts that the considered three-body 
decay processes are governed by vector-meson dominance 
(cf. fig. [J). We stress again that such a statement hinges on 
the representation used for the vector meson fields (cf. also 
the discussion in [7] ) . In our work we use tensor fields and 
our counting scheme leaves no room for contact terms at 
leading order. Other works, which use the vector represen- 
tation, do find additional contact terms, see e.g. [T^] and 
references therein. We stress again that in these works the 
experimental decay width (|ll|l is used to pin down the cou- 
pling constants of the respective model. In contrast, our 
result (fT0|) is a consequence of our effective Lagrangian 



4.1 Decay K*~ 



K+ 



We identify the three states with momenta pi, j»2, P3 with 
7r + , 7r~ , K + , respectively. The relevant matrix element to 
be used in f6} is 



Ck* + ^it+ 71— K+ 

my hp Ha 



SprriK* 

my hp m^ bA 
prriK* 



[s P {qi){ql + ™ 



L K , 






l-&f.(9i)(?i 

S K *{q{) . 



4.)] 

(14) 



The propagator of the K* meson is constructed in the 
same way as the rho- meson propagator ([9]). Also here it 
does not influence the results appreciably, if the energy 
dependent width of the K* meson is neglected. 

The Dalitz plot for this reaction is shown in fig. \5[ The 
corresponding plot with pure phase space is depicted in 
fig. O There are some deviations from pure phase space 
which are interesting to be resolved experimentally. 



The single-differential decay width 



dr K *+_ „+ 7T - k+ 



dm\ 2 



is 



plotted in fig. [7] No significant deviations from pure phase 
space remain in this integrated quantity (not shown). 
For the partial decay width one gets 



r. 



K* + ^TT+ ' 



K+ =14keV, 



(15) 



in agreement with the experimental constraints (|13p . Switch- 
ing off the effect of bA increases the width by 4 keV. The 
order of magnitude of our result (|15[) agrees with the range 
given in [13j . It clearly would be interesting to confront our 
prediction with experimental results. In high-statistics ex- 
periments like e.g. HADES at GSI and WASA at COSY 
it might be possible to observe such rare decays. 
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Fig. 5. Dalitz plot for the decay process K* — » -k it K 
(arbitrary normalization). Note ra 12 = mjj. +7r _ and m^s = 
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Fig. 6. Dalitz plot for the phase space of the decay process 

K* + ->7T + 7T- K+. 



4.2 Decay K* + -> tt° K+ it 



For the decay K* 



ir° K + 7T° we identify the three 



states 1, 2, 3 with ir°, K + , tt°, respectively. It holds 
my hp Yia 



G K . 



K* + ~*7T° K+TT° — 



16 / 3 m,K* 
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(5 



E 
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Fig. 7. Single-differential decay width for the process K* + — » 
7r + 7T~ A" 1 ". The solid line gives the full result, the dotted with 

b A = 0. 



my hp m^ bA 
2pm K . 



[S K ,(ql)-S K ,( q l)} 



(16) 



In addition, for the decay width one has to account for 
the presence of two identical states by a factor 1/2. 

The Dalitz plot for this reaction is shown in fig.[U Due 
to the destructive interference present in (fTfj]) there is a 
striking difference to pure phase space shown in fig. [3 
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, 4.3 Decay K* + -> n+ K° tt° 
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Fig. 9. Dalitz plot for the phase space of the decay process 
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Fig. 10. Single-differential decay width for the process K* — > 
7T° K + ty°. The solid line gives the full result, the dotted with 
&a = 0. For comparison also the shape for pure phase space is 
plotted as the dashed line. 



is obvious. For the partial decay width one gets the value 



Ar«+->7r° k+ tt° =67eV, 



(17) 



We identify the three states 1, 2, 3 in J6J with 7r + , K°, 
7T°, respectively and derive 
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m v hph A r „ , 2 \ I 2 , 2 
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m K . 



+ S K *{ql) {ql + m 2 K ,) + 2S p (q 2 2 ) (q 2 2 + m 2 K ,)\ 
m v h P mlb A [ ,2s, q („2\,<) m K q (t i 

(18) 



The Dalitz plot for this reaction is shown in fig. [TTJ In 
practice there is no difference to pure phase space shown 
in fig. H 
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Fig. 11. Dalitz plot for the decay process K* 
(arbitrary normalization). Note m 12 = m^+K an< ^ m 23 
m 2 K0lr0 . 



The single-differential decay width 



dr. 



K*+- 



r+ K° 7T 



dm 2 Q 



plotted in fig. [T2l For the partial decay width one gets 



r. 



K' + ^tt+ K°-n° 



28 keV. 



(19) 



whereas b A = leads to 69 eV. It has already been pointed 
out in [13] that this particular decay channel has a rather 
small width. Of course, this smallness is caused by the 
already mentioned negative interference. Hence, while a 
fully or single-differential width shows an interesting pat- 
tern, the overall missing strength of this decay channel 
probably inhibits an experimental verification. 



This result is in agreement with the experimental con- 
straints (fT3|) . We observe that for this decay the influence 
of b A is significant. For b A = the decay width comes at 
35 keV, right at the boundary of the empirical constraint 
(fT3|) . Qualitatively, also in [13] it was found that the decay 
width (|T9|) is about twice as large as the one given in (|T5|) . 
Among the three-body decay branches of the K* + me- 
son the partly charged channel 7r + K° tt° has the largest 
width. An experimental verification of our results (| 1 5|) . 
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Fig. 12. Single-differential decay width for the process K* + — » 
7r + K 7T°. The solid line gives the full result, the dotted with 
b A = 0. 



(fT7|) . Q19[) for the partial widths and maybe even for the 
differential widths shown in figs. [5J El EEQj ED] and EH 
would be highly desirable. Like in the radiative two-body 
decays studied in [3] we predict sizeable flavor breaking 
effects in the hadronic three-body decays. 



5 Summary 

The chiral Lagrangian with Goldstone bosons and light 
vector mesons has been used to calculate hadronic three- 
body decays of light vector mesons following the recently 
proposed counting rules [1] . The leading order results agree 
very well with all available data and constraints. We have 
presented predictions for differential decay widths of all 
considered processes and for the partial decay widths of 
the three-body decays of the K* meson. We suggest to 
study such reactions in high-statistics experiments like 
HADES at GSI and WASA at COSY. 
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